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Abstract 

The Statistical Model of the Early Stage, SMES, describes a transition between confined and de- 
confined phases of strongly interacting matter created in nucleus-nucleus collisions. The model was 
formulated in the late 1990s for central Pb+Pb collisions at the CERN SPS energies. It predicted 
several signals of the transition (onset of deconfinement) which were later observed by the NA49 exper¬ 
iment. The grand canonical ensemble was used to calculate entropy and strangeness production. This 
approximation is valid for reactions with mean multiplicities of particles carrying conserved charges 
being significantly larger than one. 

Recent results of NA61/SHINE on hadron production in inelastic p+p interactions suggest that 
the deconfinement may also take place in these reactions. However, in this case mean multiplicity of 
particles with non-zero strange charge is smaller than one. Thus for the modelling of p+p interactions 
the exact strangeness conservation has to be implemented in the SMES. This extension of the SMES 
is presented in the paper. 
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I. INTRODUCTION 


Strongly interacting matter at sufficiently high energy density is predicted to exist in a phase 
of quasi-free quarks and gluons, the quark gluon plasma (QGP). Relativistic nucleus-nucleus 
(A+A) collisions provide a unique opportunity to check this prediction and study properties 
of the transition to the QGP as well as the QGP itself. This is because the system created in 
A+A collisions is close to (at least local) equilibrium. The conclusion is based on the success of 
statistical and hydrodynamical models of particle production at high energies (see e.g. Ref. ED- 
Consequently, properties of the system (matter) can be characterized by its equation of state 
which should include different phases and transitions between them. It is important to note 
that nowadays there is no dynamical understanding of the observed equilibrium properties of 
particle production in A+A collisions. 

With increasing collision energy the energy density of matter created at the early stage of 
A+A collisions increases. Thus, at a sufficiently high collision energy the matter is expected 
to be created in the QGP phase. The beginning of the QGP creation with increasing collision 
energy is referred to as the onset of deconfinement. The experimental search for the onset 
of deconfinement in central Pb+Pb collisions was performed by the NA49 experiment at the 
Super Proton Synchrotron (SPS) of the European Organization for Nuclear Research (CERN) 
about 15 years ago. The study was motivated [2] by predictions of the Statistical Model of the 
Early Stage (SMES) [3] of A+A collisions. According to the model the onset of deconfinement 
in central A+A collisions should lead to rapid changes of the energy dependence of several 
hadron production properties, all located in a common energy domain. In particular, a non¬ 
monotonic dependence of the strangeness to entropy ratio as a function of the collision energy 
(the horn ) was predicted [3] as an important signal of the transition. This and other predictions 
of the SMES were confirmed by NA49 14, 5J. Moreover, following results from the Relativistic 
Heavy Ion Collider at Brookhaven National Laboratory and the Large Hadron Collider (LHC) 
at CERN agree with the NA49 results and their interpretation (see Ref. [6]). The SMES 
predictions and the experimental evidence for the onset of deconfinement are presented in 
recent reviews [7]. 

The SMES is probably the simplest model of the onset of deconfinement. This leads to 
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a number of advantages and disadvantages. In particular, the SMES is frequently criticized 
for being based on simple assumptions which can not be justified within popular dynamical 
approaches to A+A collisions. 

In this paper we concentrate on a single aspect of the SMES which concerns the finite size 
effects for strange hadron production. The SMES predictions for strangeness production were 
calculated within the grand canonical ensemble (GCE). This approximation is valid for central 
Pb+Pb collisions at the SPS energies in which mean multiplicity of particles with non-zero 
strange charge is significantly larger than one. However, this is not the case for inelastic p+p 
interactions at the SPS energies. Here the exact strangeness conservation has to be imposed 
using the canonical ensemble (CE) [13HT7]. 

Recently, the NA61/SHINE Collaboration at the CERN SPS published results on hadron 
production in p+p interactions 0 HD]. They suggest that in these reactions the strangeness 
to entropy ratio (experimentally replaced by the K + to 7r + ratio) also changes rapidly in the 
SPS energy range, see Fig. [l] However, the ratio and its energy dependence are significantly 
different from the horn measured in central Pb+Pb collisions. Can these results be explained 
by the onset of deconfinement as modelled by the SMES? The first step towards an answer 
to this question is taken in this paper by introducing the exact strangeness conservation. In 
order to allow for a direct comparison with the previously published predictions, the remaining 
SMES assumptions, parameters and notations are kept unchanged. 

The paper is organized as follows. In Sec. [TT] the GCE formulation of the SMES is briefly 


recapitulated. The exact strangeness conservation is introduced in the SMES in Sec. Ill and 


results for p+p interactions and collisions of small nuclei are presented. A summary in Sec. IV 
closes the article. 


II. THE SMES MODEL IN BRIEF 

The SMES model was formulated almost 20 years ago. Its basic assumptions, parameters 
and results are summarized in this section. Together with the notation used in the original 
paper they are here kept unchanged as much as possible in order to allow for a direct comparison 
with the previously published results. 
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Figure 1: The horn structure in the energy dependence of the K + /7r + ratio is interpreted as evidence 
for the onset of deconfinement located at low CERN SPS energies. The structure was first discovered 
by NA49 in central Pb+Pb collisions. Surprisingly its shadow is visible in inelastic p+p interactions 
as indicated by the new NA61/SHINE data. 

The SMES assumes that the matter created at the early stage of collisions has zero conserved 
charges. Consequently, its properties are entirely defined by the available energy and the 
volume in which production takes place. In central A+A collisions this volume is chosen as 
the Lorentz contracted volume occupied by the colliding nucleons (participant nucleons) from 
a single nucleus: 

v = 47rrgA p /3 
\/snn /‘ZrriN 

where m n is the nucleon mass, ^/snn is the center of mass energy of the nucleon pair, A p is the 
number of participant nucleons from a single nucleus. The ro parameter is taken to be 1.30 fm 
in order to fit the mean baryon density in the nucleus, p 0 — 0.11 fm -3 . 

Only a fraction, p, of the total energy in A+A collisions is transformed into the energy of 
new degrees of freedom created at the early stage. This is because a part of the energy is carried 
by the net baryon number. The released (inelastic) energy is expressed as 

E = rj ( y/s NN - 2 m N ) A p , (2) 
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where the parameter q is assumed to be independent of the collision energy and the system 
size. The value of q used for numerical calculations is 0.67 [3]. 

Assumptions Q and ([2]) with q — 1 correspond to the Landau hydrodynamical model HU. 
Similarly to this model, in the SMES we do not consider dynamical mechanisms leading to 
a fast thermalization of the matter. The SMES model postulates that the creation of new 
particles at the early stage of collision is a statistical process, namely, all microscopic states 
allowed by conservation laws is equally probable. 

The SMES predictions for the pion multiplicities are based on the assumption that the 
entropy generated at the early stage of collision is (approximately) conserved during the expan¬ 
sion of produced matter, ft was indeed observed that the dissipative effects estimated by the 
ratio of the shear viscosity to the entropy density are small for the strongly interacting matter, 
especially in a region of the deconfinement transition (see, e.g., [112] and references therein). 
It should be also noted that particle interactions play rather different role for the equilibrium 
properties (e.g., the equation of state) and the kinetic coefficients (e.g., the shear viscosity). 
This is clearly demonstrated by a simple example of the hard balls system H- The hard core 
particle radius r leads to small corrections to the ideal gas equation of state due to the excluded 
volume effects, but the shear viscosity as it behaves as oc r~ 2 and thus it is strongly dependent 
on r. 

The elementary particles of strong interactions are quarks and gluons. The deconhned state 
is considered to be composed of u, d and s quarks and the corresponding anti-quarks each 
with internal number of degrees of freedom equal to 6 (3 color states and 2 spin states). The 
contribution of c, b and t quarks is neglected due to their large masses. The internal number 
of degrees of freedom for gluons is 16 (8 color states and 2 spin states). The masses of gluons 
and non-strange (anti)quarks are taken to be 0. The strange (anti)quark mass is taken to be 
175 MeV [3j- The properties of equilibrated matter are characterized by an equation of state 
(EoS). For the case of quarks and gluons the bag model EoS is used [8], i.e., the ideal gas EoS 
modified by a bag constant B. This equilibrium state of quarks and gluons is called the Quark 
Gluon Plasma or Q state. 

The SMES uses an effective parametrization of the confined hadron state, denoted as W 
state. The non-strange degrees of freedom which dominate the entropy production are taken to 
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be massless bosons. Their internal number of degrees of freedom is taken to be 16 i.e., about 
3 times lower than the internal number of effective degrees of freedom in the QGP. The mass 
of strange degrees of freedom is assumed to be 500 MeV, equal to the kaon mass. The internal 
number of strange degrees of freedom is assumed to be 14. For the W-state the ideal gas EoS 
is selected. Clearly, this description of the confined state should only be treated as an effective 
parametrization. The numerical values of the parameters are fixed by fitting A+A data at the 
AGS, see for details Ref. pj. 

The model assumes that always the maximum entropy state is created at the early stage of 
A+A collisions. In the model with two different states (W and Q) the form of maximum state 
changes with the collision energy. The regions in which the equilibrium state is in the form of a 
pure W or a pure Q state, are separated by the region in which both states coexist (the mixed 
phase). The maximum entropy condition is equivalent to the assumption of the first order phase 
transitions with the Gibbs criterion for the mixed phase (see Appendix B in Ref. [3]). Namely 
at a given temperature T the system occupies a pure phase W or Q whose pressure is larger, 
the mixed phase is formed if both pressures are equal pw = Pq • The transition temperature 
between the W and Q phases is assumed to be T c = 200 MeV. 

Using the assumptions and parameters defined above predictions of the SMES can be cal¬ 
culated. The early stage energy density reads: 

_ E _ rjp 0 (y/s NN — 2m N )^/sNN 
V 2rriN 

The pressure and energy density functions in the W-phase and Q-phase are equal to: 

( 4 ) 

( 5 ) 

B, ( 6 ) 

+ B . (7) 

The strange particle contribution to thermodynamical functions (|4](7]) are taken within the 
Boltzmann approximation. This simplification is important for the CE treatment which will be 
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discussed in the next section. Note that in Ref. [3] the Fermi distribution with m*Q = 175 MeV 
was used for the strange quarks. In order to minimize differences to the previous results we 
choose here a larger value of mg = 216.5 MeV which leads to the same number of strange 
quarks at the phase transition temperature (T c = 200 MeV): 


k 2 dk exp ( — \/k 2 + ttIq/Tc ) = 


k 2 dk 


exp [\Jk 2 + m*Q /T^j + 1 


( 8 ) 


Then the bag constant B = 570 MeV/fm 3 is calculated using the Gibbs criterion of equal 
pressures: 


Pw{T c ) — Pq(T c ) . (9) 

The entropy densities in the pure phases (i =W, Q) read: 

s,(T) = pVLpm. 

In the mixed phase the W and Q phases coexist. The fraction of volume occupied 
phase is denoted as £. The energy and entropy densities in the mixed phase are 

Gni X (T C ) = £eq(T c ) + (1 — C,)£w(T c )i (11) 

Smha(T c ) = !;Sq(T c ) + (1 — £) Sw(Tc) ■ (12) 


( 10 ) 
by the Q 


The temperature T and pressure p are shown as a functions of the collision energy in Figs. [2] 
left and right , respectively. The mixed phase starts at collision energy ^sjvjv.i an d ends at 
y/SNN,2 : 

y/ snn,i = 7.42 GeV, yj Sjvjv ,2 = 10.83 GeV . (13) 

The equivalence of the Gibbs criterion and the maximum entropy condition is illustrated in 
Fig. [3]Ze/f, where the ratios R t = Si /sq are presented for i =W, mix, and Q. 

The number density of the sum of strange and anti-strange particles in the GCE can be 
calculated as 
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Figure 2: The temperature (left) and pressure (right) of the matter created at the early stage of 
A+A collisions as function of collision energy. 




Figure 3: Left: The ratio of entropy densities s^/sq with i referring to the W (solid line), Q (horizontal 
solid line), and mixed (dashed line) phases, as a function of collision energy. The full circles correspond 


to the beginning and end of the mixed phase given by Eq. (13). Right: The strangeness to entropy 
ratio n s /s as a function of collision energy. 
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In Fig. fright the strangeness to entropy ratio, n s /s, is shown as a function of the collision 
energy. 


III. PHASE TRANSITION WITH EXACT STRANGENESS CONSERVATION 


In p+p interactions at the CERN SPS energies mean multiplicity of produced strange and 
anti-strange particles is smaller than one. Thus in this case the exact strangeness conservation 
should be taken into account. In the statistical models this is done within the CE formulation. 
The CE partition function of strange particles assures an equal number of strange and anti- 
strange charges, N s = Ng, in each microscopic state of the system. For the W and Q phases it 
has a similar form and reads 

00 00 N a „Ns 1 r2n 

Z«(T,V,A)=££]v 7 ^ji(JV.-JV,)=-/ d<j> exp 

Ns=0Ns= 0 s ' s ' Jo 

where 


zie^ + e-**) = Iq(2z) , (17) 


z — Z W,Q — A V TI W q(T) . 


(18) 


The auxiliary A parameter in Eq. (18) is introduced to calculate the total strangeness density 
in the CE: 
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h Kg|T)] 
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(19) 


The ratio of Bessel functions I\ and Jo in Eq. (19) quantifies the strangeness suppression 
(relatively to the GCE yield) due to the conservation of net strangeness in each microscopic 
state of the CE. 
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In order to take into account the exact strangeness conservation for thermodynamical func¬ 
tions it is convenient to rewrite them as following: 
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where u s ^q\T, V ) is given by Eq. (19), and uw,q{T) is average energy of strange particle: 


uw.q{T) — 


J™dkk 2 (k 2 + m ,„ 0 ) 1/2 exp [- {k 2 + m f TO ) 1 / 2 /T] 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

>* 

(24) 


jy dk k 2 exp [- (k 2 + rn^q) 1 / 2 /T] 
with rri\Y taken in the W phase, and rriQ in the Q phase. The entropy density is given by 


Eq. (10). 


For A p 1 the system volume ([I|) is large, and Vn\ V Q 1. Then one finds that 


I\ [Vn^g] //q [Vn^g] —>■ 1 and, therefore, nlyj E> —> n^g. The results for the CE and GCE 


become equivalent in this thermodynamical limit, and Eqs. (20 23) coincide with Eqs 


In the mixed phase, Eq. (19) should be replaced by 


n^(T,V,0 = n‘ mQ (T) 


h[X] 

h[x] ’ 


where 


(25) 


(26) 


A' = X(T,V,0 = ZVn- q (T) + (1 - QV n‘ w (T) 

is the total GCE number of strange and anti-strange particles (both hadrons and quarks) in 
the mixed phase. This is because the CE condition of zero net strangeness in the mixed phase 
should be obeyed by the whole system and not by its phases separately. 

At each yj snn one calculates V and e according to Eqs. ([Tj) and ([3]), respectively. The CE 


Eqs. (19 23) are used to obtain results for the pure W and Q phases. In the mixed phase, the 


temperature T and the parameter £ are obtained by solving the equations: 

?£5"[T,A] + (1-Qe“ k [T,A] = e(^), 
pjf [T,X] = P r|T.A] . 


(27) 

(28) 
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where £^g and p^g are given by Eqs. (20 23) with ri^Q (25) instead of n(^g E) (19). 

The collision energy ^/Snn,i and temperature 7\, where the mixed phase starts, and yTTV/vg 
and T 2 , where the mixed phase ends, are obtained as solutions of Eqs. ( 27p8 ) for £ = 0 and 
£ = 1, respectively. One finds: 


T x = 203.4 MeV , 
T 2 = 202.9 MeV , 



y/SNN,l 

= 7.20 GeV , 

(29) 

y/SNN,2 

= 10.75 GeV . 

(30) 



Figure 4: The CE temperature (left) and pressure (right) as function of collision energy are shown 
by the solid lines. The dashed lines correspond to the GCE result presented in Fig. [2j 


The collision energy dependence of T and p obtained within the CE is shown by the solid 
lines in Fig. [4]Ze/f and right , respectively. The dashed lines in Fig. [4] correspond to the GCE 
results presented in Fig. [2] The CE and GCE curves are similar. A slightly larger value of T in 
the CE than in the GCE is needed to compensate the CE suppression of energy density. Note 
that £ as a function of yGV/v is given by Eq. (|3j) and, thus, it is independent of the system 
size. The entropy density is given by Eq. ( fT(j| in terms of p, e, and T. Therefore, the entropy 
density s is weakly affected by the exact strangeness conservation imposed in the CE. 


The Gibbs criterion (28) used in the CE is again equivalent to the maximum entropy con¬ 
dition. This is illustrated in Fig. [5]Ze/t, where the ratios R{ = (si/sq )ce calculated in the CE 


are presented for entropies sh, sq and s r , 


for A p = 1. 
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Figure 5: (a): The CE (A p = 1) ratio of entropy densities Si /sq with i referring to the W (solid line), 

Q (horizontal solid line), and mixed (dashed line) phases as a function of the collision energy. The 


full circles correspond to the beginning and end of the mixed phase given by Eq. (29) and Eq. (30), 


respectively, (b): The CE strangeness to entropy ratio as a function of the collision energy. The solid 
line corresponds to A p = 1 and the dashed line to A p > 1 which coincides with the GCE results 
presented in Fig. [3] (b). The dashed-dotted and dotted lines show the CE results for A p = 3 and 5, 
respectively. 

Figure [I] rig ht presents energy dependence of the strangeness to entropy ratio, n s /s, calcu¬ 
lated within the CE for A p = 1, 3 and 5 as well as the result for the GCE (A p 1). 

In Fig. [6] the ratio 


[n s /s\ ce 

r = 


(31) 


(»”/f>]cCE 

is shown as a function of A p at three collision energies in the vicinity of the transition region. It is 
seen that the CE suppression of the strangeness to entropy ratio depends strongly on the number 
of participants 2A p . With increasing A p the CE suppression decreases. At y / s)v)v > 10 GeV the 
suppression parameter ( [iTT] ) is close to unity already for A p > 10. The CE suppression increases 
with decreasing collision energy, when the total number of strange particles is small. This is 
indicated by the dotted line in Fig. [6] which is calculated for collisions at snn = 5 GeV. 

Finally, the strangeness to entropy ratio calculated for central Pb+Pb collisions (the GCE 
result) and inelastic p+p interactions (the CE with A p = 1 result) is plotted in Fig. [7] as 
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Figure 6: The CE strangeness to entropy ratio divided by the corresponding ratio in the GCE is 
shown as a function of A p . The double ratio (Eq. (|31[) ) is calculated at the beginning of the mixed 
phase y/sNN ~ 7.3 GeV (solid line), below the mixed phase region y / s/v r /v = 5 GeV (dotted line), and 
above the mixed phase region = 20 GeV (dashed line). 

a function of collision energy up to the LHC energies. The left plot shows the two ratios 
separately, whereas the right one presents their ratio. The energy dependence predicted by 
the SMES is only qualitatively similar to the measured one (Fig. [l|. Clearly the SMES, the 
simplest model of the onset of deconfinement, has to be significantly modified in order to reach 
a quantitative agreement with the data. 



IV. SUMMARY 

This paper introduces the exact strangeness conservation in the Statistical Model of the 
Early Stage [3j of nucleus-nucleus collisions. This allows to calculate the energy dependence of 
the strangeness to entropy ratio for collisions of protons and small nuclei at the CERN SPS 
energies. The extension of the model is motivated by the recent results of the NA61/SHINE 
experiment at the CERN SPS on hadron production in inelastic p+p interactions |10j . which 
suggest that the deconfinement may take place also in this reaction. 


13 











Figure 7: Left: Collision energy dependence of the strangeness to entropy ratio calculated within 
the SMES for central Pb+Pb collisions (the GCE result) and inelastic p+p interactions (the CE result 
for A p = 1). The ratio is plotted up to the LHC energies. Right: The ratio of strangeness to entropy 
ratios calculated for central Pb+Pb collisions and inelastic p+p interactions with the SMES as a 
function of collision energy. 

The CE treatment of the strangeness production leads to the well known effect - the total 
number of strange and anti-strange particles is reduced in comparison to that obtained within 
the GCE at the same values of volume and energy density. However, the calculations show 
only small modifications of the system temperature, pressure, and entropy density. Thus, 
the strangeness to entropy ratio is significantly reduced in small systems. The smaller the 
collision energy, the smaller is the total number of strange particles, and, thus, the stronger is 
the CE strangeness suppression. In the region of the mixed phase, y/sjy n = 7 — 11 GeV, the 
strangeness to entropy ratio in p+p interactions is found to be approximately two times smaller 
than in central Pb+Pb collisions. Note that the CE suppression becomes quite small already for 
central collisions of intermediate size nuclei and it is negligible for central Pb+Pb collisions. The 
calculated collision energy dependence of the strangeness to entropy ratio in p+p interactions 
is qualitatively similar to the one measured by the NA61/SHINE collaboration [9, 10] for the 
K + to 7r + ratio (see Fig. [Tj). However, a quantitative comparison between the model and the 
data requires further modifications of the model and thus being beyond the scope of this paper. 
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